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1.  Abstract. 


The  promise  of  new  architectures  and  more  cost-effective  miniaturization  has 
prompted  interest  in  molecular  electronics  and  molecular-based  optical  memories.  Nature 
has  already  optimized  through  serendipitous  natural  selection  some  molecules  for  such 
applications.  We  examine  here  the  use  of  the  protein  bacteriorhodopsin  in  a  branched- 
photocycle  three-dimensional  memory.  By  using  a  sequential  one-photon  process,  parallel 
read  and  write  processes  can  be  carried  out  without  disturbing  data  outside  of  the  irradiated 
volume.  This  final  report  presents  the  status  of  our  research  and  development  effort  and  an 
overview  of  what  can  be  realistically  expected  during  the  next  year  of  work  on  this  project. 

The  goals  of  this  contract  were  to  fabricate  a  prototype  protein-based  volumetric 
memory  based  on  the  branched-photocycle  architecture,  test  the  reliability  of  the  memory 
and  optimize  the  protein-based  memory  medium.  A  prototype  bench-scale  prototype  was 
developed  and  tested  to  investigate  the  viability  of  the  branched-photocycle  architecture 
based  on  bacteriorhodopsin.  The  main  difference  in  the  present  study  relative  to  previous 
prototype  studies  was  the  use  of  commercially  available  laser  diodes  to  do  both  the  paging 
and  writing  operations.  We  also  replaced  the  active  matrix  liquid  crystal  spatial  light 
modulators  with  film,  because  all  the  commercial  modulators  produced  diffraction  and 
artifacts.  A  CCD  array  detector  was  used  to  monitor  the  paged  data.  A  variety  of  data 
cuvettes  were  prepared  to  test  various  optical  geometries  and  protein  environments.  The 
protein  was  optimized  by  using  both  thermal  and  chemical  modifications,  resulting  in  a 
five-fold  improvement  in  data  write  efficiency.  Reliability  problems  remain,  but  reliability 
has  been  improved  by  a  factor  of  four  relative  to  the  earlier  Level  I  prototype.  The 
probability  of  a  read  error  is  now  0.009/bit  for  symmetric  (highly  random)  data  and 
0.011/bit  for  unsymmetric  (non-random)  data.  This  error  rate  must  be  improved 
significantly  before  this  memory  architecture  can  be  considered  to  have  commercial 
potential. 
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2.  Rationale  and  Overview 

Bacteriorhodopsin  is  the  light-transducing  protein  in  the  purple  membrane  of 
Halobacterium  salinarium,  an  archaebacterium  that  inhabks  salt  marshes  [1-5].  This 
bacterium  has  survived  five  mass  extinctions  during  a  3.5  billion  year  evolutionary  period, 
and  has  evolved  a  protein  with  high  photochemical  efficiency,  photochromic  cyclicity  and 
thermal  stability  [6-12].  A  ribbon  structure  of  the  protein  is  shown  m  Fig.  1.  In  the  native 
organism,  the  protein  converts  light  into  a  proton  gradient  that  is  used  by  the  bacterium  as  a 
chemiosmotic  source  of  energy. 

F 


Figure  1.  The  structure  of  bacteriorhodopsin  is  shown  using  ribbons  to  represent  the 
seven  alpha  helical  segments  that  span  the  membrane.  Selected  amino  acid  residues  that 
participate  in  the  proton  pumping  process  are  indicated.  The  shaded  gray  cylinder  (voluine 
5=  500  A^)  has  a  length  of  about  26  A  and  a  diameter  of  about  5  A  and  delineates  the 
approximate  volumetric  region  of  the  proton  channel. 
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Advances  in  central  processor  technology  have  altered  the  character  of  computer 
technology  from  processor-limited  towards  memory-limited  performance  [9].  The  overall 
throughput  of  many  scientific  and  image  analysis  problems  is  now  determined  primarily  by 
the  size  and  data  transfer  bandwidths  of  the  random  access  memory  rather  than  the  speed  of 
the  central  processor  unit  or  the  floating  point  hardware.  This  situation  has  generated  an 
increased  awareness  that  new  memory  architectures  providing  more  cost  effective  storage 
capacities,  local  processing  capabilities  or  data  bandwidths  should  be  investigated  [13]. 
Three-dimensional  memories  store  information  in  a  volumetric  memory  medium  and  offer 
as  much  as  a  one-thousand  fold  improvement  in  data  storage  capacity  for  a  given  enclosure 
size,  although  optical  and  reliability  considerations  tend  to  reduce  the  comparative 
advantage  factor  to  values  closer  to  300  [10,11,14,15].  The  two  most  common 
architectures  are  based  on  two-photon  [10,13-16]  or  holographic  [17-19]  methods.  The 
branched-photocycle  architecture  explored  here  is  useful  because  it  rigorously  excludes 
photochemistry  outside  of  the  doubly  irradiated  volume,  a  problem  which  must  be  solved 
before  high  density  volumetric  memories  can  reach  their  full  potential. 

The  branched-photocycle 
architecture  investigated  here  was 
developed  in  response  to  problems 
with  the  two-photon  based 
architecture.  First,  high  laser 
intensities  are  required,  which  prevent 
the  use  of  inexpensive  laser  diodes  to 
write  data.  This  observation  is  true 
even  though  bacteriorhodopsin  has 
one  of  the  largest  two-photon 
absorptivities  reported  for  a  broad¬ 
band  photochromic  material  [20-22]. 

Second,  unwanted  photochemistry 
invariably  occurs  along  the  laser  axes, 
which  requires  rather  complex 
photochemical  cleaning  operations  to 
be  carried  out  following  the  write 
process  [9,23,24].  Third,  no  protein 
state  which  provides  for  extended 
lifetime  data  storage  could  be  created 
via  a  direct  two-photon  process  [5]. 

These  three  issues  suggest  that 
commercialization  will  not  be  viable 
without  considerable  additional 
research  and  development. 

Fortunately,  during  the  course  of  these  studies  we  identified  a  branching  reaction  of  the 
protein  which  offers  a  novel  approach  to  reading  and  writing  data  within  a  three- 
dimensional  polymer  cuvette  containing  the  protein. 

The  photocycle  of  bacteriorhodopsin  at  30°C  is  shown  in  Fig.  2.  The  main 
photocycle  is  shown  in  the  left-hand  portion  of  the  figure  and  is  activated  by  the  primary 
photochemical  event,  which  converts  the  resting  state,  bR,  into  the  primary  photoproduct, 
called  K.  The  remainder  of  the  photocycle  involved  thermal  reactions  which  result  in 
reformation  of  the  resting  state,  bR,  in  about  6-10  milliseconds,  depending  upon 
temperature.  The  last  intermediate  in  the  photocycle  is  called  O,  and  this  state  can  be 
converted  by  light  to  form  P,  which  subsequently  decays  to  form  Q.  The  latter  two 
intermediates  are  the  sole  participants  in  the  branching  photocycle,  which  exists  as  an 
independent  entity  from  the  main  photocycle.  The  key  observation  here  is  that  this 
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Figure  2.  The  photocycle  of  bacteriorhodopsin. 
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branching  sequence  can  only  be  created  by  a  sequential  one-photon  process.  Furthermore, 
the  Q  state  has  a  lifetime  at  room  temperature  of  about  five  years. 

Before  we  discuss  the  details  of  the  branched-photocycle  architecture,  it  is  useful  to 
compare  the  performance  of  this  architecture  relative  to  the  two-photon  architecture  based 
on  the  use  of  the  same  protein.  A  comparison  of  the  sequential  one-photon  (branched- 
photocycle)  architecture  with  the  simultaneous  two-photon  architecture  based  on 
bacteriorhodopsin  is  shown  in  Table  I.  We  note  that  a  majority  of  the  table  entries  are 
based  on  simulations,  and  thus  we  must  delay  final  judgment  until  all  the  prototyping  work 
is  completed.  We  note,  however,  that  all  of  the  data  collected  to  date  indicate  that  our 
simulations  were  either  correct  or  conservative. 


Table  I.  Sequential  One-Photon  versus  Simultaneous  Two-Photon 

Volumetric  Memory  Architectures  Based  on  Bacteriorhodopsin(a) 


V  i  -  Archltectulreii(or«Rr^ 

TworPhoton 

Write  method 

simultaneous  2P 

Molecular  representation  of  zero  bit('^) 

bR 

bR 

Mnifiniilar  reoresentatlon  of  one  bit(^) 

M 

P,Q 

power/AT  limited 

diffraction  limited 

Maximum  write  speed  (1024x1024  pixel;  per  cube) 

Read  method 

2P  photovoltaic 

1 P  paqed  AA 

Parallel  read  capability 

power/AT  limited 

diffraction  limited 

Maximum  read  speed  (1024x1024  pixel;  per  cube) 

~1  MBvte/s 

Rit  .<;i7fi  in  cubic  microns  (1.6  cmxl  .6  cmx2  cm  cube) 

~2.3  GByte 

~1  GByte 

Conversions  per  photon 

-0.0000001 

-0.0001 

Unwanted  axial  photochemistry 

none 

Unwanted  axial  photochemistry  after  cleanina 

-2% 

n/a 

Power  required  _ 

~  1  mJ 

-0.01  mJ 

Lasers  required 

pulsed 

cw  or  laser  diode 

Siqnal-to-noise  ratio 

1.1 

1.9 

Temperature  for  >two  year  storage 

-10°C 

10-40°C 

Optimal  operatinq  temperature 

-10°C 

30°C 

Shelf  life  of  data  at  ambient  temperature 

weeks-months 

years 

(a)  romant  performance  neutral;  italics  f  acceptable  or  better; 

bold  =  acceptable  but  places  significant  limits  on  performance; 
bold  italics  =  will  prevent  commercialization  unless  improved. 

(b)  In  most  implementations,  both  the  zero  bit  and  the  one  bit  are  represented  as  linear 
combinations  of  bR  and  a  second  intermediate  (M  in  the  two-photon  memory, 
and  P:Q  in  the  sequential  one-photon  memory).  For  example,  the  branched- 
photocycle  memory  operates  with  maximal  efficiency  by  assigning  zero  to  [95% 
bR  +  5%  P+Q]  and  one  to  [85%  bR  +  15%  P+Q]. 
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3.  The  Branched-Photocycle  Architecture 

A  basic  requirement  for  the  implementation  of  all  branched  photocycle  volumetric 
memories  is  the  availability  of  a  molecule  or  light-transducing  protein  which  undergoes  a 
photocycle  upon  light  activation  and  which  has  one  or  more  thermal  intermediates  within 
the  photocycle  which  can  be  photochemically  converted  into  a  stable  species.  Orthogonal 
light  beams  are  used  to  activate  the  branched  photocycle  only  at  that  location  within  the 
volumetric  memory  cell  where  the  two  beams  cross  and  the  timing  between  the  light  beams 
is  appropriate.  We  anticipate  that  the  branched-photocycle  volumetric  architecture  can  be 
implemented  by  using  a  variety  of  different  molecules  and  light-transducing  proteins.  As 
we  demonstrate  below,  however,  bacteriorhodopsin  has  native  properties  near-optimal  for 
this  implementation.  Data  can  be  written,  read  and  stored  by  using  the  following  optical 
scheme  to  access  the  branched  intermediates,  P  and  Q. 


scheme  1 


A.  Bacteriorhodopsin  as  an  optical  AND  gate. 

Bacteriorhodopsin  (MW  =  26,000)  is  the  light  harvesting  protein  in  the  purple 
membrane  of  a  micro-organism  formally  known  as  Halobacterium  salinarium  and 
commonly  called  Halobacterium  halobium  [5,8].  The  bacterium  grows  in  salt  marshes 
where  the  concentration  of  salt  is  roughly  six  times  higher  than  that  found  in  sea  water. 
The  purple  membrane  is  grown  by  the  bacterium  when  the  concentration  of  oxygen 
becomes  too  low  to  sustain  respiration.  The  protein,  upon  the  absorption  of  light,  pumps 
a  proton  across  the  membrane  generating  a  chemical  and  osmotic  potential  that  serves  as  an 
alternative  source  of  energy.  The  fact  that  the  protein  must  survive  in  the  harsh 
environment  of  a  salt  marsh,  where  the  temperatures  can  exceed  65°C  for  extended  periods 
of  time,  requires  a  robust  protein  that  is  resistant  to  thermal  and  photochemical  damage. 
Recent  experiments  indicate  that  the  protein  can  withstand  temperatures  as  high  as  140°C 
[7],  but  under  most  conditions,  denaturation  occurs  at  80°C.  The  cyclicity  of  the  protein 
(i.e.  the  number  of  times  it  can  be  photochemical  cycled)  exceeds  10°  [5,8],  a  value 
considerably  higher  than  most  synthetic  photochromic  materials. 

The  light  absorbing  chromophore  of  bacteriorhodopsin  is  all-trans  retinal  (vitamin  A 
aldehyde).  The  chromophore  is  bound  to  the  protein  via  a  protonated  linkage  to  a  lysine 
residue  attached  to  the  protein  backbone.  The  protein  bound  chromophore  carries  a 
positive  charge,  which  interacts  elecrostatically  with  charged  amino  acids  in  the  protein 
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binding  site.  These  interactions  impart  photochemical  properties  to  the  chromophore  t^ 
differ  significantly  from  those  observed  in  the  chromophore  when  in  solution.  Upon  the 
absorption  of  light,  a  rotation  around  the  Ci3=Ci4  double  bond  occurs  to  generate  a  13-cis 
geometry  This  process  is  called  a  photoisomerization,  and  in  bacteriorhodopsin  it  is 
Complete  in  less  than  one  picosecond.  The  reason  for  this  unusual  speed  is  due  to  a 
barrierless  excited  state  potential  surface  [5].  In  this  regard,  bacteriorhodopsin  is  a 
photochemical  analog  of  a  high 

electron  mobility  transistor  'n 

(HEMT)  device  [9].  I  sesnm  I - ^  .  . 


568  nin 
paging 
laser 


„  1>  laser  Latched  i _ bR-^bR  J  out  returns  lo 

If  we  assign  the  r  and  Kl  and  1  (O^)  l  resting  state 

states  as  the  logical  state  one,  _____  gate  /  L  (returns  to  o) 

scheme  1  illustrates  the  use  of 

the  protein  as  an  optically  '  - - 

addressed  latched  AND  gate.  In  f  no  absorption 

other  words,  the  only  way  to  Latched  \ _ I  andno^change 

generate  a  logical  one  (or  in  this  I  677  nm _ I  gate  J  (stays  at  0) 

case,  to  set  bit  1 )  is  to  provide  two  d|ta  — |  ^ 

inputs  (two  photons)  at  the  correct  ■ 

wavelength  and  timing  (see  insert  553 

at  right).  This  is  the  key  to  the  i.b-,q  {  KSIUSSSsh.. 

rigorous  exclusion  of  unwanted  and  j  -  (o^i)  i  to  blue  (Q)  state 

photochemistry  and  is  responsible  ezz^nm _  gate  /  L  (converts  0  to  1) 

in  part  for  making  the  branched-  laser  I  ■**^  _ _ 

photocycle  architecture  so  I— — — — — 

attractive  in  terms  of  reliability. 

The  quantum  efficiency  of  the  .  ■ 

primary  photochemical  event  associated  with  the  page  selection  process  is  quite  high 

(0.65).  This  is  important  in  terms  of  activating  an  entire  page 

efficiency  of  the  data  write  process  is  much  lower  (0.001  -  0.02),  which  has  both 
advantages  and  disadvantages.  The  advantage  is  that  multiple  read  operations  can  be 
carried  out  without  disturbing  the  data.  The  disadvantage  is  that  longer  madiation  t™^s  ( 

4  ms)  are  required  to  generate  adequate  photochemical  conversion.  The  latter  problem, 
however,  is  circumvented  by  carrying  out  as  many  as  10^  write  operations  simultaneously 
(see  below).  Note  that  the  two  photons  are  not  absorbed  simultaneously,  and  thus  this 
memory  operates  by  a  sequential  one-photon  process,  not  a  simultaneous  two-photon 
process.  By  activating  a  page  (via  the  “page”  photon)  and  writing  the  data  (via  the  wme 
photon)  by  using  orthogonal  light  beams,  data  can  be  accessed  in  three  dimensions.  Ihe 
details  of  these  processes  are  provided  in  Section  4. 

Relatively  little  is  known  about  the  characteristics  of  the  P  and  Q  states  that 
combine  to  represent  binary  state  1 .  A  recent  study  by  Popp  et  al.  provides  evidence  to 
indicate  that  these  two  states  have  a  9-cis  conformation  [12].  Our  molecular  models  of  the 
key  states  are  shown  in  Fig.  3.  By  using  the  same  polymer  matrix  methods  as  descnbed  in 
our  two-photon  patent  [25],  volumetric  memory  media  can  be  prepared  with  optimal 
properties.  The  details  of  our  current  memory  media  are  provided  in  Section  4.D  below. 
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Figure  3.  Molecular  models  of  the  binding  site  of  bacteriorhodopsin  in  the  resting  state 
(bR  or  BIT  0),  the  P  state  (short-term  BIT  1)  and  the  Q  state  (long-term  BIT  1).  The 
location  of  the  calcium  ion  in  bR  is  based  on  two-photon  spectroscopy  [22].  The  locations 
of  the  calcium  ion  in  P  and  Q  are  not  known.  The  chromophore  is  all-tram  in  bR  and  9- 
cis  in  P  and  Q.  The  remarkable  longevity  of  the  Q  state  is  due  to  the  fact  that  the  link  to 
the  protein  is  broken,  coupled  with  the  fact  that  the  chromophore  remains  trapped  in  the 
binding  site.  The  only  quick  way  of  returning  the  Q  state  to  bR  is  to  irradiate  it  with  blue 
light,  which  isomerizes  9-cis  back  to  all-tram.  This  photoisomerization  is  followed  by 
spontaneous  reformation  of  the  chromophore-protein  link  and  regeneration  of  bR. 


4.  Implementation  of  the  Prototype 


In  this  section  we  describe  the  design  and  fabrication  of  a  bench-scale  (TyP^  I) 
prototype  and  the  methods  that  are  used  to  write,  read  and  erase  data  The  key 
modifications  in  our  prototype  relative  to  our  previous  studies  involves  the  replacement  of 
the  krypton-ion  lasers  with  commercially  available  diode  lasers.  A  schematic  diagram  of 
the  prototype  is  shown  in  Figures.  4A  and  4B.  The  write  and  read  optical  protocols  are 
summarized  in  Fig.  5  and  the  absorption  spectra  of  the  key  intermediates  m  the  mam  and 
the  branched  photocycle  are  shown  in  fig.  6. 
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Figure  4A.  Optical  diagram  of  the  protein-based  branched-photocycle  optical  memory. 
OHL,  quartz  halogen  lamp  (used  for  data  erase);  PA,  page  aperture;  DBS,  dichroic  beam 
splitter;  BEO,  beam  expanding  optics;  SLM,  spatial  light  modulator  (selects  wluch  data  bits 
within  the  page  will  be  written);  BCO,  beam  condensing  optics;  DC,  data  cuvette 
containing  the  protein  in  a  transparent  polymer  matrix;  CCD,  charge  coupled  device  (reads 
data);  DCKH,  data  cuvette  kinematic  holder;  PTC,  Peltier  temperature  controller. 

Figure  4B  (next  page):  Schematic  diagram  of  the  branched-photocycle  three-diinensiond 
memory.  The  four  operations  associated  with  the  process  of  data  storage,  retrieval  and 
erasure  are  shown.  The  paging  operation  precedes  both  the  writing  and  reading  processes; 
erasure  is  independent.  There  are  dual  paging  laser  systems  (Fig.  2),  but  we  show  only 
one  for  clarity.  The  use  of  incandescent  erasure  leads  to  multiple  pages  being  erased 
simultaneously  (see  discussion  below). 
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Figure  4B  (Caption  on  previous  page) 
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Figure  5.  A  schematic  diagram  of  the  parallel  write  and  read  sequences  with  ti^ 
progressing  from  bottom  to  top.  The  write  sequence  (left  column)  is  as  follows:  (1)  mm 
on  the  page  beam  to  activate  bR  photocycle,  (2)  after  ~2  rns,  activate  appropriate  6^m 
data  beari  to  convert  O  to  P.  The  read  sequence  (nght  column)  is  as  follows:  0)  turn  on 
568nm  page  laser  to  activate  bR  photocycle  within  paged  region  only,  (2)  rfter  ~2  ms,  reaa 
emtrpSe  as  an  image  on  the  CO)  arriy  detector.  The  basic  optical  archttectuie  is  shown 

in  insert  (a). 
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Figure  6.  Absorption  spectra  of  selected  intermediates  in  the  photocycle  of 
bacteriorhodopsin.  Gray  arrows  indicate  light  activated  processes  and  black  arrows 
indicate  thermal  reactions.  The  structures  of  the  polyene  chromophore  are  shown  as  inserts 
[N  represents  nitrogen;  X  represents  either  nitrogen  (in  P)  or  oxygen  (in  Q)]. 


A.  Data  are  written  in  parallel  via  a  branching  reaction. 

A  parallel  write  is  accomplished  by  using  the  sequential  one-photon  optical  protocol 
shown  in  Fig.  5.  The  vertical  axis  of  Fig.  5  charts  nominal  time  relative  to  the  firing  of  the 
page  addressing  (or  paging)  laser.  The  paging  beam  (k  =  630nm,  At  <  1  ms)  activates  the 
photocycle  of  bacteriorhodopsin  and  after  a  few  milliseconds  the  O  intermediate  approaches 
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maximal  concentration.  The  data  laser  and  the  SLM  are  now  activated  (X  -  690  nm,  At 
3  ms)  to  irradiate  those  volume  elements  into  which  1  bits  are  to  be  written.  This  process 
converts  O  to  P  in  these,  and  only  these,  locations  within  the  memory  cuvette.  After  mariy 
minutes  the  P  state  thermally  decays  to  form  the  Q  state  (the  P  ^  Q  decay  Xp,  s 
highly  dependent  upon  temperature  and  polymer  matrix).  The  write  process  s 
accomplished  in  ~10  ms,  the  time  it  takes  the  protein  to  complete  the  . 

represents  an  overall  write  data  throughput  of  ten  million  characters  per  second  (10  MB/s) 
for  a  1024x1024  data  array  and  two  error  correcting  bits  per  byte. 

B.  Data  are  read  in  parallel  by  using  differential  absorptivity. 

The  read  process  takes  advantage  of  the  fact  that  light  around  690  nm  is  absorbed 
by  only  two  intermediates  in  the  photocycle  of  light-adapted  bacteriorhodopsin,  the  primary 
photoproduct  K  and  the  relatively  long-lived  O  intermediate  (see  Fig  6).  A  parallel  read  is 
accornplished  by  using  a  differential  absorption  process  as  shown  in  the  right  column  of 
Fie  5  The  read  sequence  starts  out  in  a  fashion  identical  to  that  of  the  write  process  by 
activating  the  568nm  paging  beam.  After  two  billiseconds,  data  laser  is  activated  and  the 
data  read  (DRS)  shutter  is  opened  for  1ms,  but  the  SLM  is  left  off  allowing  only  0.1%  of 
the  total  laser  power  through.  A  CCD  array  (clocked  to  clear  all  charges  prior  to  opening 
DTS  and  DRS)  images  the  light  passing  through  the  data  cuvette.  Those  elements  m 
binary  state  1  (P  or  Q)  do  not  absorb  the  690nm  light,  but  those  volumetric  elements  that 
started  out  in  the  binary  0  state  (bR)  absorb  the  690nm  light,  because  these  elements  have 
cycled  into  the  O  state.  Noting  that  all  of  the  volumetric  elements  outside  of  the  paged  area 
are  restricted  to  the  bR,  P  or  Q  states,  the  only  significmt  absorption  of  the  beam  is 
associated  with  O  states  within  the  paged  region.  The  CCD  detector  pay  therefore 
observes  the  differential  absorptivity  of  the  paged  region,  and  the  papd  region  alone.  This 
selectivity  is  the  key  to  the  read  operation,  and  it  allows  a  reasonable  signal-to-noise  ra  m 
even  with  thick  (1-1.6  cm)  memory  media  containing  >103  pages.  Because  the 
absorptivity  of  the  O  state  within  the  paged  region  is  more  than  1000  tmes  larger  than  the 
absorptivity  of  the  remaining  volume  elements  combined,  a  very  weak  beam  can  be  used  to 
generate  a  large  differential  signal.  The  read  process  is  complete  in  ~10  ms  which  gives  a 
rate  of  10  MB/s.  Each  read  operation  must  be  monitored  for  each  page,  and  a  refresh 
operation  performed  after  -1000  reads.  While  data  refresh  slows  the  memory  slightly, 
page  caching  can  minimize  the  impact. 

C.  Data  can  be  erased  by  page  or  globally. 

In  our  previous  (Level  I)  prototype,  we  used  a  Krypton-ion  laser  operating  at  413 
nm  to  erase  single  pages  at  a  time.  Our  current  prototype  (Level  HA)  is  designed  to  test  the 
potential  of  using  commercially  available  diode  lasers,  but  here  we  haw  a  problem.  At 
present  there  are  no  readily  available  blue  diode  lasers,  although  a  number  of  companies 
^e  working  on  developing  these  lasers.  In  the  interim,  we  are  using  a  filtered  quartz- 
halogen  lamp  (QHL  in  Fig.  4A)  to  provide  the  blue  light  for  erasure.  By  choosing  an 
incandescent  source,  any  hope  of  single  page  erasure  is  lost,  and  now  we  must  be  satisfied 
with  multiple  page  erasure.  For  the  time  being,  we  accept  this  constrained,  anticipating  tha 
blue  diode  lasers  will  be  available  in  the  near  ftiture. 

The  optimal  wavelength  for  erasing  data  is  -410  nm,  because  this  wavelength  will 
photoconvert  both  P  and  Q  back  to  bR  (scheme  1  &  Fig.  6).  Alternatively  one  can  c  ear 
an  entire  data  cuvette  by  using  incoherent  light  in  the  360  -  450  nm  range.  The  latter  option 
may  prove  useful  for  some  commercial  implementations.  A  broad-band  data  clear 
operation  must  be  carried  out  in  the  absence  of  red  light,  however,  because  blue  light  can 
activate  the  photocycle  and  produce  the  O  state.  As  long  as  the  erasing  light  is  n^ 
absorbed  appreciably  by  the  O  state,  however,  a  sufficient  population  of  the  P  and  Q 
molecules  (bit  1)  will  be  reset  to  bR  (bit  0). 
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As  an  aside,  it  is  interesting  to  consider  whether  a  clear  operation  should  be 
implemented  at  all.  The  cost  of  a  data  cuvette  is  relatively  small.  As  noted  below,  we  use 
inexpensive  plastic  cuvettes  and  the  amount  of  protein  that  is  required  to  make  a  single 
cuvette  capable  of  storing  hundreds  of  megabytes  is  on  the  order  of  milligrams.  Because 
the  protein  can  be  produced  in  mass  quantities  by  using  fermentation  technology  not  much 
different  from  making  beer,  some  investigators  have  proposed  that  the  memory  be 
implemented  as  a  WORM  device  (write  once-read  many  times).  An  inexpensive  bulk 
eraser  device  using  incoherent  blue  light  could  be  provided  for  those  wishing  to  re-use  data 
cuvettes. 

D.  Preparation  of  the  data  cuvette. 

Bacteriorhodopsin  (BR)  is  isolated  from  the  bacterium  Halobacterium  halobium 
(also  called  Halobacterium  salinarium)  in  the  form  of  purple  membrane  sheets,  which 
consist  of  the  protein  suspended  in  a  lipid  membrane.  Our  current  strain  of  H.  halobium  is 
a  variant  of  S9-P,  an  enhanced  overproducer  of  bacteriorhodopsin,  which  is  isolated  via 
well-defined  procedures  described  in  the 
literature  [22,26].  Purple  membrane  sheets 
are  suspended  in  distilled  water  and  stored 
at  4°C  until  ready  for  use.  Prior  to  use,  the 
concentration  of  BR  is  adjusted  to  10-20 
mg/ml. 

The  cubes  used  in  the  branched 
photocycle  volumetric  memory  are 
produced  by  suspending  purple  membrane 
in  a  polymer  matrix.  The  matrix  must  be 
sufficiently  rigid  to  encapsulate  the 
membrane  and  trap  it  in  place,  while 
providing  enough  moisture  to  ensure  that 
bacteriorhodopsin  retains  its  full  activity.  It 
also  must  provide  a  high  degree  of  optical 
clarity  and  homogeneity  for  proper 
functioning  of  the  optical  memory.  The 
polymer  matrix  must  be  rigid  enough  to 
prevent  protein  migration,  but  flexible 
enough  to  permit  the  protein  to  undergo 
modest  changes  in  volume  that  take  place 
during  the  photocycle.  Although  we  are 
still  looking  for  improved  polymer 
matrices,  to  date  poly(acrylamide)  has 
shown  the  most  promise. 


Fabrication  of  bR-based  optical 
cuvettes  (Fig.  7)  is  accomplished  by  in  situ 
polymerization  of  the  acrylamide-bR 
solution.  Concentrated  purple  membrane  is 
added  drop-wise  to  a  small  volume  of  10  % 

(wt/vol)  acrylamide  (~2  ml)  until  the  optical 
density  at  570  nm  reaches  2-3  absorbance  units.  The  resulting  solution  is  then  filtered  (5 
p,,  syringe-type)  and  sonicated  for  1-2  minutes  (prolonged  sonication  can  result  in  vesicle 
formation,  inadequate  sonication  generates  a  polymer  with  too  much  scattering).  At  this 
point  the  solution  is  degassed  for  several  minutes,  because  dissolved  oxygen  can  inhibit  the 
acrylamide  polymerization  process.  Fresh  ammonium  persulfate  solution  [(NH4)2S208, 
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Figure  7.  Schematic  diagram  of  the 
protein  data  cuvette  and  holder. 
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1.5%  aqueous,  w/v)]  is  added  to  initiate  and  catalyze  polymerization,  followed  by  a  small 
volume  of  TEMED  (N,  N,  N',  N'-Tetramethyl-ethylenediamine)  to  act  as  an  accelerator  of 
the  polymerization  process.  The  volumes  of  these  two  reagents  are  typica  ly  on  the  order 
of  75  uL  and  10  iiL,  respectively,  per  mL  of  10%  acrylamide.  The  solution  must  be 
quickly  mixed  and  transferred  to  a  standard  fluorescence  cuvette  before  polymerization  is 
fully  accelerated;  failure  to  do  so  will  result  in  optical  distortions  within  the  matrix  (floxvs 
and  eddies  due  to  differing  rates  of  polymerization).  The  of  f  largely 

controlled  by  the  concentrations  of  ammonium  persulfate  and  TEMED.  The  voluiries  of 
these  reagents  can  therefore  be  adjusted  to  minimize  formation  of  optical  defects  in  the 
cube.  The  total  volume  of  acrylamide/bR  solution  needed  is  approximately  1.5  ml. 

The  acrylamide  solution  is  prepared  to  be  10%  T  and  5%  C,  where  T  represents  the 
total  monomeric  concentration  of  acrylamide,  and  C  represents  the  coricentration  of  cross- 
linking  agent  (bisacrylamide).  Addition  of  a  cross-linking  agent  enables  the  formation  of 
an  actual  gel,  as  opposed  to  a  viscous  solution  of  long  polymer  chains  which  would  be 
incapable  of  encapsulating  the  purple  membrane  fragments.  The  relative  proportions  of  T 
and  C  are  crucial  to  the  ultimate  optical  quality  of  the  gel,  especially  with  respect  to  light 
scattering.  The  concentration  of  cross-linking  agent  controls  the  effective  pore  size,  pd  is 
generally  within  the  range  of  2.5  -  5%  (pore  size  of  approximately  20  nm).  Addition  of 
too  much  bisacrylamide  (relative  to  acrylamide)  leads  to  enlarged  pore  diameters  (several 
hundred  nm)  and  a  consequent  increase  in  light  scattering  [27]. 


After  the  cube  has  gelled,  a  second  acrylamide  layer  is  added.  The  ^rpose  of  this 
layer  is  two-fold:  primarily,  it  aids  in  the  reduction  of  light  scattering  and  reflecticm  which 
might  interfere  with  the  reading  and  writing  processes,  and  secondly  it  acts  as  a  filler  or 
the  rest  of  the  cuvette,  which  is  then  sealed  to  ensure  that  it  is  air-tight.  The  protocol  tor 


fabrication  of  this  layer  is  as  described 
above,  except  that  enough  bR  is  added 
so  that  the  solution  appears  black: 
deeply  colored  laser  dyes  were  originally 
attempted  for  this  purpose,  but 
preliminary  studies  indicated  that  they 
interfered  with  the  polymerization 
process.  The  cuvette  is  filled  to  about 
two  thirds  of  its  total  volume,  and  is 
made  air-tight  by  the  use  of  either 
silicone  sealants  or  a  UV-curable  optical 
adhesive.  A  custom-made  base  for  the 
cube  facilitates  permanent  sealing  of  the 
cuvette,  as  well  as  temperature 
regulation  via  pel  tier  coolers.  The  cube 
base  is  illustrated  in  Fig.  7.  The  cuvette 
is  inserted  into  the  holder,  which  is 
equipped  with  ports  for  insertion  of 
thermocouples  which  monitor  the 
temperature  and  indirectly  control  the 
peltier  coolers. 


Figure  8.  The  effect  of  temperature  on 
the  magnitude  and  the  decay  time 
of  the  differential  read  signal. 


E.  The  effect  of  temperature. 

Our  studies  indicate  that  the  signal-to-noise  ratio  of  the  read  process  is  highly 
dependent  upon  temperature.  Furthermore,  the  probability  of  writing  a  one  bit  is  enhanced 
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by  increasing  the  temperature.  Both  factors  are  associated  with  the  effect  of  temperature  on 
the  formation  of  the  O  state.  The  key  experimental  results  are  shown  in  Fig.  8.  The 
downward  peak  is  associated  with  the  formation  of  the  O  state,  and  the  consequential 
absorption  of  the  676.5  nm  beam  by  this  state.  (We  note  that  our  present  prototype  uses  a 
690nm  diode  laser,  but  the  data  shown  in  Fig.  8  are  relevant  to  any  read  operation  in  the 
650  -  700  nm  region.)  The  rise  time  is  associated  with  the  decay  of  the  O  state,  but  the  fact 
that  it  "over  shoots"  the  zero  level  is  due  to  capacitive  effects  in  the  detector  array  (we  are 
working  to  eliminate  this  problem).  Note  that  the  intensity  of  the  negative  signal  as  well  as 
the  decay  of  the  signal  increases  with  increasing  temperature.  This  observation  can  be 
explained  by  reference  to  the  detailed  kinetic  analyses  of  the  photocycle  carried  out  by  Varo 
and  Lanyi  [28,29].  The  signal  to  noise  ratio  of  the  read  process  is  proportional  to  the 
magnitude  of  the  negative  signal  divided  by  the  magnitude  of  the  signal  fluctuation  in  the 
time  segment  from  -2  to  0  ms.  The  worst-case  S/N  ratios  for  the  read  process  as  a  function 
of  temperature  are  as  follows:  1.5  (10°C),  3  (20°C),  8  (30°C)  and  11  (40°C).  Other 
sources  of  error  become  dominant  in  producing  read  errors  at  temperatures  above  30°C, 
and  hence  we  have  adopted  this  temperature  as  the  nominal  operating  temperature  of  the 
data  cuvette.  As  noted  previous,  the  temperature  of  the  cuvette  is  controlled  by  using 
thermoelectric  heaters  (see  Fig.  7). 

Temperature  also  affects  the  timing  of  the  write  and  read  data  beams.  For  optimal 
writing,  the  on  and  off  times  are  given  by  the  following  equations. 

At  (data  beam  on)  =  21.8  ms  Tc(°C)-9-7  (10°C<Tc<45°C)  (1) 

At  (data  beam  off)  =  61 . 1  ms  Tc(°C)-0-7  (10°C  <  Tc  <  45°C)  (2) 

where  At  is  the  time  in  milliseconds  relative  to  address  beam  activation  and  Tc  is  the 
temperature  of  the  memory  cuvette  in  degrees  Centigrade.  The  above  relationships  are 
approximate  and  are  limited  to  pH=7  solution  or  well  humidified  polymer  environments, 
fte  data  beam  on  time  is  equal  to  the  time  at  which  the  signal  curve  if  Fig.  8  reaches  half  of 
its  minimum  value.  The  data  beam  off  time  is  equal  to  the  time  at  which  the  signal  returns 
to  its  half-minimum  value.  These  times  are  appropriate  for  both  the  read  and  write  data 
beam  timings,  but  as  noted  above,  the  read  process  is  carried  out  by  using  a  very  low  light 
level.  Examination  of  Eqs.  1  and  2  indicates  that  the  memory  can  operate  more  rapidly  at 
higher  temperatures,  but  temperatures  above  40°C  should  be  avoided  as  elevated 
temperatures  decrease  the  cyclicity  of  the  protein. 


F.  Reliability 

As  noted  in  our  previous  progress  report  (RL-TR-96-274),  the  key  goal  of  the 
present  contract  was  to  improve  upon  the  reliability  of  the  read  operation.  This  has  been 
accomplished,  as  shown  in  Fig.  9.  Note  that  the  overlap  region  has  become  much  smaller 
relative  to  the  inte^al  under  the  correct  main  bit  0  and  bit  1  curves.  If  the  assignment  level 
is  set  at  the  arrow  in  each  example,  the  error  associated  with  the  read  process  is  minimized. 
In  the  following  discussion,  we  provide  a  detailed  discussion  of  the  methods  and 
procedures  of  calculating  read  error  rates  from  the  read  histograms. 


Figure  9.  (next  page)  Histograms  of  bit  intensities  during  read  operations  for  the  previous 
Level  I  prototype  (top  graph)  and  present  Level  HA  prototype  (middle  and  bottom  graphs). 
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Figure  9  (Caption  on  previous  page) 
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Calculation  of  Reliability  from  Read  Histograms.  The  improvement  in  the 
reliability  of  our  Level  IIA  prototype  is  sufficient  to  warrant  a  more  detailed  discussion  of 
this  issue,  and  how  one  can  use  the  data  presented  in  Fig.  9  to  calculate  the  probability  of  a 
read  error.  We  note  at  the  outset  that  we  place  all  of  the  error  into  the  read  portion  of  the 
process,  although  the  error  could  be  due  entirely  to  an  imperfect  write  operation.  For  the 
purposes  of  analysis,  however,  the  error  does  not  show  up  until  one  tries  to  carry  out  a 
read  operation,  and  thus  it  is  more  useful  to  analyze  error  entirely  in  terms  of  "read  error". 

The  process  of  assigning  the  state  of  a  bit  during  the  read  process  is  based  on  the 
measurement  of  the  intensity  of  light  passing  through  that  microscopic  portion  of  a  page 
upon  which  the  bit  has  been  written.  In  this  discussion,  we  assign  Bit  0  to  the  bR 
molecule  and  Bit  1  to  the  P/Q  states.  Thus,  Bit  0  is  assigned  based  on  the  preferential 
absorptivity  of  the  O  state  (for  details,  see  Figs.  5  and  6  and  related  discussion). 
Naturally,  we  could  reverse  the  bit  assignments  without  affecting  the  reliability,  and  there 
may  be  some  instances  when  such  a  reversal  would  be  useful.  At  this  stage,  we  also 
neglect  reliability  enhancement  associated  with  the  addition  of  error  correcting  bits.  This 
will  be  the  subject  of  the  final  paragraph  of  this  section. 

Whether  a  given  location  is  assigned  to  Bit  0  or  Bit  1  is  determined  by  a  threshold 
level  called  the  "bit  assignment  crossover  point",  and  represented  by  the  symbol,  The 
optimal  value  is  determined  by  carrying  out  a  reliability  analysis  on  the  assignment  of  both 
bits,  under  the  assumption  that  one  seeks  a  point  which  minimizes  the  error  in  arbitrary 
assignments.  In  other  words,  the  error  is  minimized  for  simultaneous  assignments  of  both 
bits  0  and  1  within  a  given  page.  To  understand  the  process,  and  the  assignment  of  error, 
in  more  detail,  we  provide  the  analyses  presented  in  Figures  10  and  11. 

Our  first  example  assumes  identical  Gaussian  distributions  for  bit  response 
functions  based  on  a  normalized  response  function,  from  0.0  to  1.0  (see  Fig.  10,  top). 
The  error  in  the  read  process  is  numerically  assigned  by  defining  the  following  integrals: 

Ao  =  integral  under  Bit  0  read  histogram  =  J  Bo[^]d^  (1) 

0 

1 

Aox  =  Bit  0  error  integral  =  I  Bo[^]d^ 

1 

Ai  =  integral  under  Bit  1  read  histogram  =  J  B 

0 

Aix  =  Bit  1  error  integral  =  J  B  i  [^]d^ 

0*^ 


Aoi  =  integral  under  Bit  0  and  Bit  1  overlap  regions  = 

0 


(2) 

l[^]d^  (3) 

(4) 

1 

J  Boi[|ld4  (5) 
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Figure  10.  Error  analysis  for  two  identical  overlaping  Gaussian  distributions. 
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Figure  11.  Error  analysis  for  two  different  overlaping  Gaussian  distributions. 


In  the  example  shown  in  Fig.  10,  we  have  chosen  arbitrarily  a  full-width  at  half¬ 
maximum  that  produces  an  overlap  in  the  read  histograms  (distributions).  This  overlap 
region  generates  error,  because  a  read  operation  that  returns  a  value  for  the  response 
function  within  this  region  is  ambiguous.  The  ambiguity  can  be  minimized  by  judicious 
choice  of  the  threshold  value,  C,  but  never  eliminated  completely. 

The  error  in  reading  Bit  0,  (as  a  function  of  0  is: 

^o(C)  =  in  Bit  0  read  process  =  Aq^/Aq  (6) 

with  a  similar  equation  defined  for  Bit  1: 

§l(Q  =  error  in  Bit  1  read  process  =  Aix/Ai  (7) 

The  optimal  value  for  the  bit  assignment  threshold  (or  crossover  point)  is  given  by  that 
value  of  C  for  which  io(0  =  ^l(0-  This  approach  is  shown  graphically  in  Fig.  10, 
bottom,  and  in  the  case  of  the  two  symmetric  Gaussians,  yields  ^  =  0.5  and  a  read  error  of 
eo(C=0.5)  =  <Si(C=0-5)  =  0.00927.  This  means  that  there  is  a  0.9%  chance  of  an  incorrect 
bit  assignment  during  a  single  read  process.  If  one  reads  a  page  of  data  comprising 
1024x1024  =  1,048,576  bits,  9720  bits  will  be  missassigned.  This  is  clearly 
unacceptable,  but  as  noted  below,  the  use  of  error  correcting  bits  modifies  the  situation 
significantly. 

The  reader  may  wonder  how  much  math  is  actually  involved  in  analyzing  the 
equality,  So(0  =  as  a  function  of  C-  Unfortunately,  such  calculations  normally 

require  numerical  methods,  which  are  quite  cumbersome  and  a  powerful  workstation  is 
required  if  a  large  error  histogram  data  base  is  being  used.  If  the  histograms  can  be  fit  to 
Gaussians,  the  following  pseudo-closed-form  expressions  can  be  used  to  calculate  the 
integrals: 


Let  B[^]  =  exp(  ’ 

where  cl  is  a  constant  (-4*ln(2))  yields  normalized  distributions),  xq  is  the  center  of  the 
Gaussian  distribution  and  Ax  is  the  full-width  at  half-maximum.  The  general  inte^al,  from 
a  to  b  within  the  normalized  response  function  can  be  calculated  using  the  expression, 

b 

Ax  7cU2  erf(ci^^^  Ax'^  (b  -  xq))  Ax  erf(ci  Ax~^  (q  —  xq)) 

=  2ciU2  2  Cl  1/2 

is  the  error  function.  This  function  is  normally  provided  in  newer 
implementations  of  FORTRAN  and  C  and  is  a  built-in  function  in  Mathematica  and  Maple. 

If  this  function  is  not  available,  one  is  better  off  doing  nunierical  integration  using  the  raw 
histogram  file.  It  should  further  be  noted  that  normalization  of  the  histograrn  data  is  not 
required  for  the  above  functions  to  work.  Normalization  is  automatically  provided  by  Eqs. 

6  and  7.  It  is  however  important  that  the  Bit  0  and  Bit  1  histograms  be  scaled  the  same 
amount,  if  scaling  is  being  carried  out.  We  use  normalized  data  solely  for  the  purposes  of 
plotting. 


J  Bo[^]d^ 

a 

where  "erf" 
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Given  the  computational  complexity  of  using  trial  and  error  methods  to  solve  ^o(C) 
=  <Si(Q,  we  mention  a  simple  method  of  estimating  the  minimum  error  without  assigning 
This  method  only  requires  integration  of  the  overlap  region  and  ratioing  of  that  integral 
to  the  integrals  under  the  0  and  1  bit  histograms; 


Aqi 

Aq  +  Ai 


(8) 


If  the  data  are  perfectly  symmetrical,  this  estimate  is  identical  to  the  calculated  value  based 
on  solution  of  the  §o(C)  ~  equation.  For  example,  Ae  data  in  Fig.  10  yield  0.00927 
based  on  simultaneous  minimization  of  the  <$o(C)  ®i(0  Eq.  8  gives  the  same 

result. 


Aoi  0.00395 

Ao  +  Ai  “  0.2129  +  0.2129 


=  0.00927 


Although  the  above  (Soi-based  method  works  well  for  error  calculations  involving 
symmetric  distributions  (identical  but  shifted  bit  0  and  1  histograms),  this  method  is  not 
rigorous  when  uneven  distributions  are  analyzed.  As  shovm  in  Fig.  1 L  an  unsymmetrical 
overlap  region  requires  the  explicit  evaluation  of  the  ©o(0  ^^d  ©i(0  functions  to 
rigorously  determine  the  error.  In  particular,  Eq.  8  yields  0.0462  for  the  estimated  error 
(see  data  in  top  graph  of  Fig.  1 1)  compared  to  a  more  rigorous  value  of  0.0497  (bottom  of 
Fig.  11). 


Reliability  analyses  on  the  previous  and  current  prototypes  are  summarized  in  Figs. 
12-14.  The  results  are  tabulated  below: 


Bit  Distribution 

Minimum  Error 

Fig.  No. 

Level  I 

roughly  even 

0.524 

0.044 

12 

Level  ELA 

roughly  even 

0.768 

0.009 

13 

Level  HA 

more  O's  than  I's 

0.744 

0.011 

14 

We  note  that  while  the  error  has  improved  by  a  factor  of  at  least  four,  the  read  error  is  still 
relatively  high  and  would  preclude  commercialization,  if  further  improvement  were  not 
possible.  For  reference,  we  note  that  the  simple  ©oi  "based  error  calculations  yield  values 
of  0.041, 0.009,  and  0.007  for  the  three  cases  shown  in  Figs.  12,  13  and  14,  respectively. 
Thus,  the  simple  error  calculation  breaks  downs  significantly  when  uneven  bit  distributions 
are  encountered.  Of  greater  importance  to  the  implementation  of  the  memory  is  the 
observation  that  the  optimal  bit  assignment  threshold,  shifts  in  value  when  the  ratio  of  0 
to  1  bits  changes.  This  observation  may  require  that  the  data  for  each  page  be  manipulated 
to  maintain  a  more  even  distribution,  which  is  certainly  possible,  but  rather  undesirable. 
An  alternative  is  to  use  gray  scale  recording  which  enhances  both  data  storage  density  and 
helps  to  even  out  the  distribution.  Much  remains  to  be  done  before  gray-scale  recording  is 
possible,  however. 

In  closing,  we  note  that  the  addition  of  error  correction  bits  dramatically  improves 
reliability.  For  example,  by  using  two  error-correcting  bits  for  each  8-bit  byte  of  binary 
data,  one  achieves  single  bit  error  correction  and  double  bit  error  detection.  By  using  a  10- 
bit  byte,  therefore,  we  can  decrease  the  probability  of  a  catastrophic  read  error  for  a  byte  of 
information  from  ~8%  (0.01*8)  to  less  than  0.6%.  Clearly,  error  correction  helps,  but  not 
enough  to  transform  our  present  level  of  reliability  into  commercial  viability.  More  work  is 
needed  (see  Section  7  below). 
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Normalized  Frequency  of  Observation 


Probability  of  Bit  Assignment  Error  Normalized  Frequency  of  Observation 
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Figure  13. 

Read  error  analysis  for  the  level  IIA  prototype  developed  as 
part  of  this  research  contract  for  data  sets  involving 
roughly  equal  numbers  of  0  and  1  bits. 
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Figure  14. 

Read  error  analysis  for  the  level  IIA  prototype  developed  as 
part  of  this  research  contract  for  data  sets  involving 
significantly  more  0  bits  than  1  bits. 
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5.  Protein  Modifications  to  Improve  Memory  Performance. 

The  storage  of  data  using  P  and  Q  states  is  linked  directly  to  O  state  yield,  and  we 
have  been  studying  both  chemical  and  genetic  methods  of  improving  this  variable.  Most  of 
this  work  was  carried  out  during  the  previous  contract,  and  we  summarize  the  current  state 
of  knowledge  here.  One  promising  approach  is  site  directed  mutagenesis,  and  in  particular, 
the  replacement  of  leucine  93  with  tyrosine  or  the  replacement  of  glutamic  acid  204  with  a 
non  carboxylate  amino  acid  [30-32].  These  single  amino  acid  substitutions  greatly  increase 
the  yield  of  the  O  state,  but  they  also  increases  the  lifetime  of  the  O  state,  which  slows 
down  the  cycle  time  of  the  memory.  Other  mutants  are  under  study.  This  section  will 
concentration  on  a  discussion  of  organic  cation  protein  analogs,  a  few  of  which  yield  a 
five-fold  improvement  in  protein  performance. 

It  has  long  been  known  that  the  functioning  of  light-adapted  bacteriorhodopsin 
(^ax  ®  570  nm)  as  a  proton  pump  requires  the  presence  of  Ca2+  and  Mg2+  ions  which 
bind  at  specific  locations  inside  the  protein  and  on  the  surface  [22].  Approximately  3-5 
divalent  cations  are  bound  to  the  protein,  although  the  mechanism  of  the  binding  remains  a 
subject  of  debate  [33].  Removal  of  these  cations  produces  what  is  called  the  blue  membrane 
(^max  ~  605  nm),  a  protein  with  a  truncated  photocycle  incapable  of  vectorial  proton 
pumping.  The  spectroscopic  and  photochemical  properties  of  the  protein  can  be  restored 
fully  by  adding  either  calcium  or  magnesium  to  an  aqueous  solution  of  the  blue  membrane 
[22,33].  The  second  high-affinity  binding  site  is  primarily  responsible  for  the  blue  to 
purple  transition  [34].  Recent  two-photon  experiments  indicate  that  this  binding  site, 
hereafter  referred  to  as  the  chromophore-adjacent  cation  binding  site,  involves  the 
participation  of  ASPgs  and  ASP212  along  with  additional  stabilization  from  TYR57  and 
TYRi85  [22].  This  model  of  the  binding  site  is  shown  in  the  bottom  diagram  of  Fig.  3. 

We  have  observed  that  analog  proteins  with  nearly  identical  spectroscopic 
properties  and  similar  photochemical  properties  can  be  generated  by  adding  large  organic 
monovalent  and  divalent  cations  to  the  blue  membrane  [35].  A  selection  of  some  of  the 
organic  cations  that  we  have  studied  with  respect  to  protein  performance  related  to  the 
three-dimensional  memory  are  shown  in  Fig.  15.  The  spectra  of  the  organic-cation  bR 
analogs  are  similar  to  the  native  protein  with  absorption  maxima  within  a  few  nanometers 
of  the  native  protein.  In  order  to  examine  systematic  variations,  the  main  absorption 
bands  of  the  protein  spectra  were  fit  by  using  a  log-normal  distribution  function.  This 
procedure  fits  the  band  to  a  skewed  G^ssian  as  a  function  of  the  absorption  maximum 
(Xmax).  full-width-at-half-maximum  (Av)  and  skewness  (a).  The  results  of  the  least- 
squares  fits  are  shown  in  Fig.  15  for  spectra  recorded  in  polyacrylamide  gel  (the  polymer 
used  to  make  the  data  cuvettes).  As  we  demonstrate  below,  the  modest  shift  in  absorption 
spectra  do  not  fully  reveal  the  large  changes  that  were  observed  in  the  O  state  kinetics. 
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Figure  15.  Structure  of  the  quaternary  ammonium  cations  and  the  spectroscopic 
properties  of  the  light-adapted  analog  proteins  generated  via  addition  of  the  organic  cations 
to  the  blue  membrane  of  bacteriorhodopsin. 


A.  Preparation  of  the  organic  cations. 

The  dibromoalkanes  and  trialkyl  amines  were  purchased  from  Aldrich  (cat.  nos. 
24,065-6,  12,590-3  and  14,080-5)  and  Fluka  (cat.  nos.  92260,  90342  and  93240) 
respectively  and  were  used  as  received.  Acrylamide,  N,N'-methylene-bis-acrylamide, 
azolectin,  tris(hydroxymethyl)aminomethane  (Tris)  and  EDTA  were  purchased  from  Sigma 
(cat.  nos.  A-8887,  M-7279,  P-7443  ,T-1503  and  EDS  respectively)  and  were  used 
without  further  purification.  Ammonium  persulfate  and  CaCl2.2H20,  purchased  from 
Fisher  (cat.  nos.  A682-500  and  C-79  respectively)  were  of  pure  grade.  lodomethane  and 
tetramethylethylenediamine  were  purchased  from  Aldrich  (cat.  nos.  28,956-6  and  41,101- 
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9  respectively)  and  were  used  as  received.  Distilled  and  deionized  water  (Barnstead  E- 
pure)  was  used  throughout  the  experiments. 

The  divalent  ‘bolaform’  cations  [CnR'6  where  n  is  the  number  of  carbon  atoms  in 
the  alkyl  chain  attached  between  the  two  quaternary  ammonium  groups  and  R'  references 
the  six  identical  hydrocarbon  groups  attached  to  the  two  nitrogen  atoms,  R'  =  CH3  (Me), 
C2H5  (Et)  or  C3H7  (Pr)]  were  synthesized  from  the  appropriate  dibromoalkane  by  using 
the  following  reaction  to  form  the  desired  n-alkane-a,co-bis(trialkylammonium)dibromide; 


MeOH 


Br(CH2)„Br  +  2  R3N  - >  BrR3N(CH2)„NR3Br  (scheme  2) 

reflux  ~3hrs. 


where  n  =  3  or  4  and  R  =  CH3  (Me),  C2H5  (Et)  or  C3H7  (Pr).  The  details  of  the 
preparation  follows  a  standard  method  [36].  A  two-fold  excess  of  the  trialkylamine  of 
choice  is  added  to  a  methanolic  solution  of  the  dibromoalkane  and  refluxed  for  about  three 
hours.  The  product,  which  is  precipitated  upon  cooling  the  solution,  is  filtered  off  and 
washed  three  times  with  cold  methanol.  The  supernatant  is  partially  evaporated  by  using  a 
rotary  evaporator  and  a  few  drops  of  ether  is  added  before  the  solution  is  chilled  to  form  a 
second  batch  of  the  product.  The  crude  product  is  recrystallized  three  times  from 
methanol  to  produce  the  pure,  white-crystalline  powder.  We  discovered,  however,  that  the 
above  procedure  is  not  applicable  to  the  preparation  of  C2Me6.  The  NMR  spectrum 
revealed  displacement  of  only  one  of  the  bromides  of  1 ,2-dibromoethane  by  trimethyl 
ammonium.  Thus,  C2Me6  was  prepared  by  the  following  procedure.  A  slight  excess  of 
methanolic  iodomethane  is  added  dropwise  into  a  methanolic  solution  of 
tetramethylethylenediamine.  The  mixture  is  left  to  stand  for  ~1  hour.  The  crude  product 
which  precipitates  is  filtered  off,  washed  three  times  with  cold  methanol  and  recrystallized 
three  times  from  methanol  to  yield  the  desired  white-crystalline  product. 

B.  Photophysical  properties  of  the  organic  cation  analog  proteins. 

Although  the  organic  cations  have  an  impact  on  a  number  of  the  photophysical 
properties  of  the  protein,  we  will  concentrate  on  an  analysis  of  only  two  properties  relevant 
to  this  contract.  The  first  is  cyclicity,  or  the  number  of  times  the  protein  can  be  cycled  from 
the  resting  state  (bR)  to  the  branched  states  (P  &  Q),  and  then  back  again  to  bR.  The 
number  for  the  native  protein  is  about  10^,  which  means  the  native  protein  under  normal 
applications,  will  function  for  about  five  years  of  normal  use.  All  of  the  organic  cation 
analogs  appear  to  have  a  comparable  cyclicity  with  the  exception  of  the  large  monovalent 
cation  analogs  such  as  BU4N.  The  latter  has  a  cyclicity  of  about  10^,  which  is 
unacceptable.  We  believe  the  latter  cation  is  too  large  to  fit  inside  the  cation  binding  site 
without  disturbing  the  location  of  the  membrane  spanning  helices.  In  the  process,  the 
protein  is  destabilized,  which  limits  the  number  of  times  the  protein  will  undergo  a 
photocycle  prior  to  denaturation. 

The  key  property  provided  by  some  of  the  organic  cations  involves  the  lifetime  and 
the  yield  of  the  O  state.  The  experimental  results  are  shown  in  Fig.  16. 
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Figure  16.  The  influence  of  organic  cation  size  and  structure  on  the  formation,  decay  and 
yield  of  the  O  state.  Some  organic  cations  enhance  significantly  the  yield  and  the  lifetime 
of  the  O  state.  The  efficiency  of  photogenerating  the  P  state  from  the  O  state  (i.e.,  writing 
bit  1)  is  proportional  to  the  integral  of  O  concentration  with  respect  to  time,  which  is 
graphed  as  a  function  of  the  cation  molecular  weight  in  (b). 


The  results  shown  in  Fig.  16  indicate  that  three  organic  cations  generate  a  dramatic 
improvement  in  the  O  state  yield  (monovalent  species  Me4N  and  Et4N,  and  divalent 
species  CaMeg;  see  Fig.  15  for  the  structures).  In  particular  Et4N  produces  a  5.4-fold 
increase  in  O  state  yield,  which  translates  to  a  potential  improvement  of  P  &  Q  formation 
from  3%  to  16%.  Our  simulations  indicate  that  optimal  performance  is  obtained  at 
approximately  10%  conversion,  so  other  organic  cations  will  also  be  investigated  to  find 
the  best  compromise  involving  yield,  stability  and  synthetic  convenience. 
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The  key  advantage  of  using  this  chemical  approach,  as  opposed  to  genetic 
engineering,  is  that  the  extent  of  the  O  state  population  and  lifetime  enhancement  can  be 
tailored  accurately  to  achieve  a  specific  goal.  This  observation  follows  from  the  fact  that 
there  is  a  very  large  number  of  possible  organic  cations,  and  by  varying  either  structure 
and/or  molecular  weight,  systematic  control  can  be  achieved. 


Figure  17.  Molecular  models  of  the  bR,  M  and  O  states  of  bacteriorhodopsin  showing 
the  chromophore,  selected  amino  acid  residues  and  the  near-chromophore  caleium  binding 
site.  Direction  towards  the  cytoplasmic  and  extracellular  surfaces  and  the  proton  pumping 
vector  are  indicated.  Microwave  absorptivity  experiments  suggest  that  in  native  bR,  the 
calcium  enters  the  proton  channel  and  provides  an  electrostatic  gate  preventing  back  transfer 
of  the  proton  during  the  latter  stages  of  the  photocycle  [37].  The  arrow  in  the  O  state 
diagram  indicates  the  migration  of  the  cation  back  to  its  native  position  during  the  O  —>  bR 
reaction. 


We  close  this  section  with  an  hypothesis  on  how  the  organic  cations  are  affecting 
the  yield  of  the  O  state.  Molecular  modeling  suggests  that  the  O  bR  reaction  requires  a 
migration  of  the  cation  (calcium  in  the  native  protein)  as  shown  in  Fig.  17.  Thus,  as  the 
molecular  weight  of  the  cation  increases,  the  decay  of  the  O  state  is  slowed.  Thus,  an 
increase  in  molecular  weight  of  the  cation  will  increase  the  yield  of  the  O  state  for  kinetic 
reasons.  However,  at  some  point  the  cation  becomes  too  large  to  occupy  the  native  cation 
binding  site,  and  thus  an  alternative  site  is  found  that  is  less  restrictive.  We  believe  this  site 
is  closer  to  the  cation  binding  site  shown  for  bR  (see  discussion  of  the  two  models  of  the 
O  state  cation  binding  sites  in  Ref.  [37]).  Thus,  there  is  a  maximum  in  the  yield  versus 
molecular  weight  graph  (Fig.  15b). 
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6.  Comments  and  Conclusions 


The  prototype  shown  in  Fig.  4  was  redesigned  to  permit  the  use  of  commercially 
available  laser  diodes  and  CCD  array  detectors.  The  goal  was  to  make  a  small-scale  level  11 
prototype  to  explore  more  fully  the  commercial  viability  of  our  architecture,  examine  better 
methods  of  paging  and  test  new  protein  analogs  developed  in  the  latter  stages  of  this 
contract.  The  key  observation  was  a  factor  of  four  improvement  in  the  reliability  of  the 
read  operation  as  discussed  in  Section  5E.  We  have  not  as  yet  found  a  spatial  light 
modulator  that  has  sufficient  reliability  and  dynamic  range  to  function  adequately  in  our 
design.  We  have  not  had  the  financial  resources  to  purchase  and  try  the  ferroelectric 
modulators  that  are  currently  being  developed  at  the  Univ.  of  Colorado.  We  hope  to  do  so 
in  the  future. 

The  protein  bacteriorhodopsin  has  many  characteristics  that  are  near-optimal  for  use 
in  the  branched-photocycle  architecture.  However,  it  is  unrealistic  to  expect  that  nature  has 
optimized  the  protein  for  this  particular  application,  and  we  will  continue  our  investigation 
of  chemically  and  genetically  modified  bacteriorhodopsin  proteins  to  identify  candidates 
that  improve  performance.  While  the  organic  cation  analog  proteins  show  great  promise, 
we  have  not  fully  evaluated  the  potential  of  site-directed  mutagenesis  to  enhance  the 
branching  photochemistry.  Through  a  collaboration  with  Prof.  Richard  Needleman,  we 
are  initiating  a  study  of  selected  mutants  with  extended  O  state  lifetimes  and  yields  to 
determine  the  impact  of  these  mutations  on  the  performance  of  the  protein  in  the  branched- 
photocycle  architecture.  Recent  studies  indicate  that  the  O  state  lifetime  and  population  can 
be  increased  significantly  by  using  Glu-204  and  Leu-93  mutants  [30-32].  When  combined 
with  careful  pH  control  [38]  and  the  use  of  organic  cations  (see  above),  we  anticipate  that 
we  can  enhance  the  amount  of  P  and  Q  that  are  formed  by  an  order  of  magnitude.  This 
will  improve  reliability  and  data  write  speeds. 

Finally,  the  homogeneity  of  the  memory  medium  requires  iniprovement.  While 
polyacrylamide  provides  excellent  protein  stability,  the  small  cavities  that  are  created 
following  polymerization  provide  scattering  centers  that  limit  resolution  and  the  quality  of 
the  read  signal  image.  DuPont  Corporation  has  offered  to  help  us  identify  new  polymers 
that  will  provide  adequate  stabilization  of  the  protein  while  improving  homogeneity  of  the 
data  cuvettes.  Hybrid  polymers  involving  various  combinations  of  polyvinyl  alcohol  and 
gelatin  will  be  investigated.  We  are  working  with  DuPont  scientists  currently.  We  have 
also  established  a  collaboration  with  one  of  Canada's  leading  young  polymer  chemists. 
Prof.  Wayne  Wang  at  Carleton  University  (Ottawa).  He  has  already  demonstrated  that  a 
co-polymer  decreases  scattering.  We  are  in  the  process  of  setting  up  a  formal  collaboration 
with  Prof.  Wang  to  explore  the  use  of  co-polymers  in  stabilizing  the  protein  in  3D  memory 
cuvettes. 
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